Collective excitations of a trapped boson-fermion mixture across demixing 
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We calculate the spectrum of low-lying collective excitations in a mesoscopic cloud formed by 
a Bose-Einstein condensate and a spin-polarized Fermi gas as a function of the boson-fermion re- 
pulsions. The cloud is under isotropic harmonic confinement and its dynamics is treated in the 
collisional regime by using the equations of generalized hydrodynamics with inclusion of surface 
effects. For large numbers of bosons we find that, as the cloud moves towards spatial separation 
(demixing) with increasing boson-fermion coupling, the frequencies of a set of collective modes show 
a softening followed by a sharp upturn. This behavior permits a clear identification of the quan- 
tum phase transition. We propose a physical interpretation for the dynamical transition point in a 
confined mixture, leading to a simple analytical expression for its location. 

PACS numbers: 03.75.Kk, 03.75.Ss, 67.60.-g 



I. INTRODUCTION 

After the achievement of Bose-Einstein condensation 
in alkali-atom gases, advanced techniques are being de- 
veloped to cool gases of fermionic isotopes inside mag- 
netic traps. Since the Pauli principle forbids s-wave 
collisions between spin-polarized fermions, to reach the 
degeneracy regime in the fermionic case one must re- 
sort to collisions against a distinguishable species, either 
bosonic or fermionic Q] . Boson-fermion mixtures are cur- 
rently being produced and studied in several experiments 

Starting from the work on strongly interacting 3 He- 
4 Hc liquids, cold mixtures have played an important role 
in the development and testing of the theory of quan- 
tum phase transitions. In this context trapped mixtures 
of atomic gases offer a unique opportunity due to their 
high diluteness. Mixtures of a Bose-Einstein condensate 
and a degenerate Fermi gas are predicted to have a rich 
phase diagram. In the case of attractive boson-fermion 
interactions, where the boson-fermion overlap is largest, 
the boson-induced fermion-fermion attraction may lead 
to the formation of a superfluid state In the case 
of repulsive interactions, on the other hand, the system 
is expected to undergo spatial separation when the re- 
pulsions overcome the kinetic energy . The conditions 
for demixing have been derived for a homogeneous mix- 
ture [lfj and also in the experimentally relevant case 
of a mixture under harmonic confinement [ill Il2t Il3| . 
The static equilibrium properties of such mixtures across 
phase separation and the topology of the particle density 
profiles in the demixed state have been studied [l4|- Of 
course, the transition to the demixed state in a meso- 
scopic cloud under confinement is spread out as the over- 
lap energy between its two components reaches a maxi- 
mum and then gradually decreases on further increase of 
the boson-fermion coupling. 

In the mixed state the dynamical properties of a boson- 
fermion mixture have been investigated both for a ho- 
mogeneous system and in a mesoscopic cloud under 
external harmonic confinement The purpose of the 



present paper is to follow the dynamics of a harmon- 
ically confined cloud with increasing boson-fermion re- 
pulsion and to look for a signature of the transition to 
spatial separation in the spectrum of collective modes. 
We focus here on the study of the dynamics in the col- 
lisional regime, which at high temperature has already 
been reached in some experiments 7| and can be attained 
at low temperature in the presence of impurities • We 
solve the hydrodynamic equations beyond the Thomas- 
Fermi approximation, including surface effects which are 
crucial for a proper description of particle density fluctu- 
ations as the cloud approaches demixing. 

The paper is organized as follows. In Sec.|n]we intro- 
duce the specific system that we study and the hydro- 
dynamic equations which we use throughout the paper. 
Section lTTTl discusses the equilibrium density profiles that 
are needed to evaluate the collective modes of the mix- 
ture in Sec. II VI Finally, Sec. ]V\ presents a summary of 
our results and an outlook towards future developments. 



II. THEORETICAL MODEL 

We consider a dilute fluid composed by two species 
of alkali atoms, one fermionic and the other bosonic 
in a Bose-Einstein condensed state, confined inside a 
spherical trap at zero temperature. The interactions 
between the bosons and between bosons and fermions 
are described by contact potentials and are parametrized 
by the coupling constants gBB = inh' Obb/ttib and 
9b F = 2irh 2 aBF/'ni r in terms of the s-wave scattering 
lengths cibb and clef and of the masses tub and tuf 
of each species, with m r = (1/itib + being 
the reduced mass. The fermions are spin-polarized and 
are taken as noninteracting, since collisions in the s-wave 
channel arc forbidden by the Pauli principle. In the fol- 
lowing we have chosen qbb > and gsF > 0, as for the 
6 Li- 7 Li mixture studied in the experiments of Schreck et 
al. 0. 

We describe the dynamics of the system by starting 
from the equations of generalized hydrodynamics [l8j for 
the particle densities p a (r,t) and the current densities 
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j<r(r, i), with a = B,F. These equations read 



dt 



+ V-j CT = 



and 



ot 



(1) 



(2) 



Here, V a are the effective mean-field potentials and W 
are the kinetic stress tensors. In the dilute regime we may 
adopt the Hartree-Fock approximation for the effective 
potentials, 
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where V£ xt = m a ui^r 2 /2 are the (isotropic) external 
trapping potentials. 

The above equations can be closed in the collisional 
regime, where we assume a local dependence of the stress 
tensors on the particle densities. In the dilute limit the 
Thomas-Fermi approximation yields the fermionic stress 
tensor as the local-density form of the tensor for the ideal 
Fermi gas. However, we have added to this Thomas- 
Fermi form a surface contribution in the form derived 
by von Weizsacker [19| , in order to avoid spurious diver- 
gences in the density fluctuations at the classical radius 
of the cloud. Thus, the form for lif^ reads 
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where A = ?i 2 (67r 2 ) 2 / 3 /2mp. This choice is in agree- 
ment with the general structure of the fermionic stress 
tensor under harmonic confinement as demonstrated in 
Ref. [2(j. In the same approximation the bosonic stress 
tensor has only the surface contribution 
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as can also be obtained from the Gross-Pitaevskii equa- 
tion |2l|. In Eqs. (J5J and JBJ we have left out velocity- 
dependent terms which do not enter linear dynamics. 

Using Eqs. ©-© for the effective potentials and the 
kinetic stress tensors we can rewrite the equations for the 
current densities as 



dt 



= P o (P CT - W CT - 9BF Vp 5 ) , 



(7) 



where a denotes the component different from a and for 
convenience we have introduced the forces 
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and 
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In these equations we shall set p a (r, t) — p a (r) + 5p a (r , t) 
and proceed first to discuss the equilibrium profiles /c CT (r). 



III. EQUILIBRIUM PROFILES AND SPATIAL 
SEPARATION 

The particle density profiles at equilibrium are ob- 
tained by imposing the steady-state condition dj a /dt = 
in Eq. J7J. This ensures consistence between static and 
dynamical solutions as well as fulfilment of the general- 
ized Kohn theorem 22]. The above equilibrium condition 
is in fact equivalent to a minimization of the mean-field 
energy functional 
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where the quantum pressure or surface energy terms read 
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6 



2m B 



6 mj? 



'Pf\ 



(11) 



(12) 



As was found in previous studies 0, 0, ITU IT^ . on 
increasing the boson-fermion repulsion the mixture un- 
dergoes spatial separation. In a finite cloud the tran- 
sition is smooth and can be described by following the 
behavior of the boson-fermion interaction energy Ei nt = 
9bf J d 3 r pbPf, which for a given couplin g st rength is 
determined by the overlap of the two species [3, Q . The 
onset of the transition is signalled first by a decrease of 
Ei n t as a function of the boson-fermion coupling (partial 
separation), until the value E int becomes negligible (full 
demixing). The maximum of the interaction energy as 
a function of gs f has been estimated for harmonic con- 
finement at zero temperature within the Thomas-Fermi 
approximation to lie at 
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and d = (Ti/mB^B) 1 ^ 2 is the bosonic harmonic-oscillator 
length. The point of full demixing as given by the 
Thomas-Fermi approximation occurs instead at 

»S? = (S)" 2 (15) 

where k F = (48N F ) 1/6 /d and a = [3 1 / 3 /(2w) 2 / 3 ](m F + 
itib) 2 /(imFrriB)- As we shall see below, the analysis of 
the mode frequencies as functions of the boson-fermion 
coupling strength yields a dynamical condition for spatial 
separation which is intermediate between the two static 
conditions laid out in Eqs. H13|l and 1)150. 



IV. COLLECTIVE EXCITATION SPECTRUM 

We proceed to evaluate the dynamical behavior of the 
cloud as it undergoes spatial separation. Under the as- 
sumption of a weak external drive we can neglect an- 
harmonic contributions. The spectrum of collective ex- 
citations is then obtained by linearizing Eq. (J7J around 
the equilibrium state. By taking the Fourier transform 
with respect to the time variable we thus obtain coupled 
eigenvalue equations for the density fluctuations 5p a of 
each species, 

m a n 2 5p a = V ■ {p a <5F CT ) -QbfV- {p a ^5pa) ■ (16) 

The expressions for the linearized forces ^F CT are given in 
Appendix IaI 

We have solved Eq. (|16fl in the case of a spherically 
symmetric confinement, setting the values of the trapping 
frequencies at ujp = u>b = with lu = 2ir x 1000 s~ 4 
and the values of the boson-boson scattering length at 
a bb = 0.27 nm, which correspond to the 6 Li- 7 Li mix- 
ture in the Paris experiment 0, H3| ■ We vary the mu- 
tual scattering length cibf and the number of particles 
of each species in order to explore the transition from the 
mixed state to the fully separated state. In particular, 
the choice Np = 10 4 and Nb = 2.4 x 10 7 leads to bosonic 
and fermionic clouds having approximately the same size 
at zero coupling. 

The numerical procedure that we have used can be 
summarized as follows: (i) we look for a spherically sym- 
metric steady-state solution of Eqs. 10-10 by a steepest- 
descent method 0H1|; we decompose the density 
fluctuations into components of definite angular momen- 
tum, i.e. we factorize the amplitude of the fluctuations as 
<5p CT (r) = 6p l a (r)Yi m (f) with Yi m the spherical harmonic 
functions; and (Hi) we set up an eigenvalue problem for 
each / by discretizing Eqs. I|16|l and solve it by means of 
standard routines from the LAPACK library |26| . 

The results for the frequencies of the monopolc (I = 0) 
modes are given in Fig. ^ as functions of the boson- 
fermion coupling for Np — 10 4 and various values of 
Nb- A sampling of the eigenvectors is reported in Fig. 
|5J We have labelled the modes as "fermionic" (dots) or 



"bosonic" (circles) according to the nature of their eigen- 
value in the limit of vanishing gBF- Of course, at finite 
values oi gBF the modes are coupled and the labels are 
just conventional, but can still be assigned by looking at 
the nodes of each density fluctuation: fermionic (bosonic) 
modes keep a constant number of nodes in Spp (8pb) 
with increasing gsF- 

A common feature of our results in Fig. ^ is a non- 
monotonic behavior of the frequency of the fermionic 
modes as the cloud evolves from the mixed regime to 
the fully separated one. With increasing gsF a fermionic 
mode acquires the character of an out-of-phase fluctua- 
tion of the two components (see Fig. |2J) and we observe 
a softening of its frequency. This we interpret as a signal 
of the approaching spatial-separation transition: with in- 
creasing boson-fermion repulsion the out-of-phase oscil- 
lation requires less and less energy, until at the transition 
the mixture takes as its equilibrium configuration a state 
which corresponds to a "frozen" out-of-phase oscillation. 
The nonzero value of the lowest mode frequency at the 
transition point appears to be due to the presence of the 
confinement, since in a linearized theory it should tend 
to zero in the proper thermodynamic limit. 

The transition point as dynamically determined by the 
sharp upturn in the fermionic mode frequencies does not 
agree with the static criterion ljT3|l or (|T5|l . but corre- 
sponds to an intermediate point where the equilibrium 
density of the fermionic cloud vanishes at the center of 
the trap (see the dotted lines in Fig. [21 . In the limit of 
large number of bosons this point can be analytically es- 
timated within the Thomas-Fermi approximation to be 



CIBB TUB + mp [LB 

where pp and ps are the chemical potentials for fermions 
and bosons, respectively. In a local-density picture this 
corresponds to a minimum in the velocity of the fermionic 
sound wave. The locations of the transition point accord- 
ing to the criterion in Eqs. i|13|) and l|17f) are shown by 
the arrows in Fig. ^ 

On further increasing the boson-fermion coupling in 
the spatially separated regime we observe that the fre- 
quency of the fermionic modes continues to increase. 
This can be understood by means of a simple model of 
sound-wave propagation inside a uniform shell of given 
thickness r\. The velocity c s of the wave is related to 
the compressibility k of the fermionic gas according to 
c s = (rappp k) -1 / 2 and by imposing rigid boundary con- 
ditions at the edges of the shell we can obtain its fre- 
quency s pect rum. For I = the lowest sound mode is 
given by [27[ 

fisw = — (18) 
V 

and the sound velocity can be estimated from the ideal- 



4 



gas expression of the pressure (p = | Ap 5 J 3 ) as 



Since we are assuming a constant-density shell for each 
value of a^f , to actually compare the predictions of 
Eq. I|18|) with the numerical results in Fig. ^ we need 
a suitable choice of the effective fermionic density pp. In 
Fig. we show the frequency of the lowest sound mode 
calculated from Eq. i|18|) as a function of cibf using two 
choices for the effective density, one corresponding to the 
maximum value of the fermionic equilibrium profile and 
the other evaluated for a uniform spherical shell with the 
same number of fermions and the radii taken from the 
Thomas-Fermi equilibrium profile. In view of the rough 
approximations that we are making, the agreement in 
Fig. |3 between the numerical calculation of the lowest 
monopole frequency and this simple model is quite sat- 
isfactory. Such a degree of quantitative agreement is not 
found for the high-frequency modes, since the density 
fluctuation tails become more marked and are more sen- 
sitive to the inhomogeneity of the fermionic shell. 

Returning to Fig. 2] it also shows that the frequen- 
cies of the bosonic monopole modes are essentially unaf- 
fected by the boson-fermion interaction. This is due to 
the fact that with our choice of system parameters the 
bosonic cloud has a considerably higher density than the 
fermionic one, leading to a weakly coupled dynamics. It 
is worth recalling that this is at present a relevant exper- 
imental situation, with the fermionic cloud containing a 
denser Bose-Einstein condensate. 

We have also examined the case of comparable num- 
bers of bosons and fermions (see Fig. 0J. In this case 
the frequency of the low-lying bosonic monopole modes 
increases slowly with increasing boson-fermion coupling 
while the frequency of the fermionic ones is quite un- 
changed except for level crossings (left panel in Fig. . 
The dipolar (1 = 1) modes show similar features (right 
panel). The behavior found for the higher modes, e.g. 
those with frequencies larger than approximately 7 loq , is 
instead similar to that shown in Fig. ^ an d is therefore 
not reported in Fig. 0] These features can be understood 
by considering the relative size of the two clouds: the hole 
produced in the fermionic cloud by the bosons is in this 
case so small that the low-lying fermionic fluctuations 
cannot sense it. 

Finally, we have solved the hydrodynamic equations 
for the I = 1 dipolar oscillations with a choice of particle 
numbers analogous to that made for the monopole modes 
in Fig. ^ As is shown in Fig. [3J in this case the lowest 
fcrmionic-modc frequency initially decreases monotoni- 
cally and tends to a constant value for large values of 
gBF- hi addition, we find the Kohn mode at the fre- 
quency of the trap for any value of the boson-fermion 
coupling. The higher modes exhibit instead the same 
behavior as the monopole modes in Fig. ^ 



V. SUMMARY AND CONCLUDING REMARKS 



In summary, in this paper we have studied how the 
mutual repulsive interactions affect the spectrum of col- 
lective excitations in a trapped boson-fermion mixture 
in the collisional regime as the mixture undergoes spa- 
tial separation. For this purpose we have derived and 
solved the equations of generalized hydrodynamics be- 
yond the Thomas- Fermi approximation, by including sur- 
face density-gradient terms in the form first proposed by 
von Wcizsacker. When the two component clouds have 
similar sizes (implying Nb ^> Np), we have found that 
the frequencies of the fermionic I — modes decrease as 
the boson-fermion scattering length is increased and the 
mixture approaches demixing. This frequency softening 
is directly related to the change in shape of the equilib- 
rium density profiles as the repulsive interactions become 
stronger and reflects the tendency of the two components 
to spatially separate. At the point where the fermions 
are expelled from the center of the trap, although the 
two clouds still partially overlap, the fermionic mode 
frequencies start to grow in a fashion which essentially 
agrees with a simple model of sound wave propagation 
inside a fermionic spherical shell. A similar trend is 
also found in the frequencies of the fermionic I = 1 bulk 
modes, while the bosonic dipolar surface mode displays 
the Kohn-theorem behavior and is therefore unaffected 
by the interactions. All these features are very differ- 
ent from those that we observe for Nb — Np, when 
the bosonic modes are the most sensitive to the boson- 
fermion coupling. 

We have thus found a clear dynamical signature of the 
onset of spatial demixing in the spectrum of collective 
modes in the case Nb ^> Np. This is expected to be 
helpful since the formation of a symmetric "egg" config- 
uration in the demixed cloud is not easily detected from 
an analysis of column density profiles [l3| | . The dynam- 
ical condition for demixing is given by the point where 
the topology of the fermionic equilibrium density profile 
changes as the fermions start to arrange themselves in 
a shell around the bosons. This point does not coincide 
with the points of partial or full demixing as obtained 
from the static study of the boson-fermion energy func- 
tional. However, this is easily understood if one considers 
that in a finite system the transition is smooth and can be 
characterized by several different conditions, which will 
coincide only in the thermodynamic limit. 

The present analysis can be extended to study the 
transition towards "exotic" configurations of the demixed 
cloud, as predicted by the static study of the energy func- 
tional |14|. A calculation of the spectrum of collective 
excitations of a boson-fermion mixture across spatial sep- 
aration in the collisionless regime is in progress and will 
be reported elsewhere. 
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APPENDIX A: FULL EXPRESSIONS FOR 
FORCE FLUCTUATIONS 

As was mentioned in the text, we decompose the parti- 
cle density fluctuations in their angular-momentum com- 
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In these equations a prime means a derivative with re- 
spect to r. 



[1] B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 
[2] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, 

T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 

Lett. 87, 080403 (2001). 
[3] J. Goldwin, S. B. Papp, B. DeMarco, and D. S. Jin, Phys. 

Rev. A 65, 021402 (2002). 
[4] Z. Hadzibabic, C. A. Stan, K. Dicckmann, S. Gupta, 

M. W. Zwierlein, A. Gorlitz, and W. Ketterle, Phys. Rev. 

Lett. 88, 160401 (2002). 
[5] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 

Rev. Lett. 89, 150403 (2002). 
[6] G. Modugno, G. Roati, F. Riboli, F. Ferlaino, R. J. 

Brecha, and M. Inguscio, Science 297, 2240 (2002). 
[7] F. Ferlaino, R. J. Brecha, P. Hannaford, F. Ri- 
boli, G. Roati, G. Modugno, and M. Inguscio, cond- 

mat/0211051. 

[8] M. J. Bijlsma, B. A. Heringa, and H. T. C. Stoof, Phys. 
Rev. A 61, 053601 (2000); D. V. Efremov and L. Viverit, 
Phys. Rev. B 65, 134519 (2002). 
[9] K. M0lmer, Phys. Rev. Lett. 80, 1804 (1998); N. Nygaard 
and K. M0lmer, Phys. Rev. A 59, 2974 (1999). 
[10] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A 

61, 053605 (2000). 
[11] A. Minguzzi and M. P. Tosi, Phys. Lett. A 268, 142 



(2000). 

[12] T. Miyakawa, K. Oda, T. Suzuki, and H. Yabu, J. Phys. 
Soc. Jpn. 69, 2779 (2000). 

[13] Z. Akdeniz, P. Vignolo, A. Minguzzi, and M. P. Tosi, J. 
Phys. B 35, L105 (2002). 

[14] Z. Akdeniz, A. Minguzzi, P. Vignolo, and M. P. Tosi, 
Phys. Rev. A 66, 013620 (2002). 

[15] S. K. Yip, Phys. Rev. A 64, 023609 (2001). 

[16] T. Miyakawa, T. Suzuki, and H. Yabu, Phys. Rev. 
A 62, 063613 (2000); T. Miyakawa, T. Suzuki, and 
H. Yabu, ibid. 64, 033611 (2001); P. Capuzzi and 
E. S. Hernandez, ibid. 64, 043607 (2001); T. Sogo, 
T. Miyakawa, T. Suzuki, and H. Yabu, ibid. 66, 013618 
(2002); P. Capuzzi and E. S. Hernandez, J. Low. Temp. 
Phys. 126, 425 (2002). 

[17] M. Amoruso, I. Meccoli, A. Minguzzi, and M. P. Tosi, 
Eur. Phys. J. D 8, 361 (2000). 

[18] N. H. March and M. P. Tosi, Ann. Phys. (NY) 81, 414 
(1973). 

[19] C. F. von Weizsacker, Z. Phys. 96, 431 (1935). 

[20] A. Minguzzi, P. Vignolo, and M. P. Tosi, Phys. Rev. A 

63, 063604 (2001). 
[21] A. Minguzzi, M. L. Chiofalo, and M. P. Tosi, Phys. Lett. 

A 236, 237 (1997). 



6 



[22] E. Zaremba and H. C. Tso, Phys. Rev. B 49, 8147 (1994); 

J. F. Dobson, Phys. Rev. Lett. 73, 2244 (1994). 
[23] Regarding the choice of the traps we have set the bosonic 

and fermionic frequencies equal, taking the geometrical 

average of the experimental values in the three Cartesian 

directions. 

[24] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes in Fortran 77: The Art of 
Scientific Computing (Cambridge University Press, Cam- 
bridge, 1992). 

[25] F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 



(1996). 

[26] E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Demmel, 
J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, 
A. McKenney, et al., Lapack User's Guide (1999), URL 
http: //www .netlib . org/lapack 

[27] This is readily obtained by imposing rigid boundary con- 
ditions on the lowest I = solution for sound-wave prop- 
agation inside a uniform spherical shell with hard walls. 
The wave is described by a spherical Bessel function jo 
of appropriate argument. 



0.4 0.8 1.2 1.6 

a RF [nm] 




FIG. 1: Frequencies of I = collective modes (in units of the trap frequency u>o) as functions of the boson-fermion scattering 
length a,BF (in nm). The panels from top to bottom correspond to Nb = 2.4 x 10 7 , 10 6 , and 10 5 bosons. Each frame displays 
the low-lying hydrodynamic modes (fermionic - dots and bosonic - open circles) from the numerical solution of Eq. 1161 with 
Nf — 10 4 fermions and the other parameters as specified in the text. In each panel the small arrow indicates the onset of 
partial demixing from a static criterion (Eq. I|l.'-^ 1 and the large arrow indicates the point of vanishing fermionic density at the 
center of the trap (Eq. 1171 1. The lines are a guide to the eye. 
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FIG. 2: Monopolar density fluctuation profiles (arbitrary units) as functions of the radial coordinate (in units of the harmonic 
oscillator length d) for various values of the boson-fermion coupling as indicated in the panels, with Nb = 2.4 x 10 7 and 
Nf — 10 4 . Continuous and dashed lines display the fermionic and bosonic profiles, respectively. The scale of each fluctuation 
has been changed to fit in the graph (F/B indicates the ratio between the maximum value attained by the fermionic and 
bosonic fluctuation profiles). The dotted lines show the equilibrium density profiles, again rescaled to fit in the graph. 




FIG. 3: Frequency of the lowest fermionic monopole mode after demixing, as a function of clbf (in nm) for Nb = 2.4 x 10 7 
and N F = f0 4 . The circles are from the numerical solution of Eq. 1161 , while the lines show the predictions of Eq. 118H taking 
as effective density the maximum value of the equilibrium density (full line) or the density of a uniform shell with Nf = 10 4 
fermions inside the Thomas- Fermi radii (dashed line). 
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FIG. 4: Frequencies of the I = (left panel) and I = 1 (right panel) collective excitations (in units of the trap frequency a>o) 
as functions of abf (in nm) for the case Nb = Nf = 10 5 . The symbols and the system parameters are as in Fig. 



a RP [nm] 



FIG. 5: Frequencies of the / = 1 modes (in units of the trap frequency u>o) as functions of clef (in nm) for Nb 
and Nf = 10 4 . The symbols and the system parameters are as in Fig.0 



= 2.4 x 10 7 



